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We propose a broadly applicable formalism for the description of coarse grained en- 
I tropy production in quantum mechanical processes. Our formalism is based on the Husimi 

transform of the quantum state, which encodes the notion that information about any quan- 
tum state is limited by the experimental resolution. We show in two analytically tractable 
cases (the decay of an unstable vacuum state and reheating after cosmic inflation) that the 
growth rate of the Wehrl entropy associated with the Husimi function approaches the classi- 
cal Kolmogorov- Sinai entropy. We also discuss various possible applications of our formalism, 
^ l' including the production of entropy in the early stages of a relativistic heavy ion collision. 

^ '. 

, ■ §1. Introduction 

I The production of entropy at the quantum level, i. e. the loss of information 

about the state of the system under consideration, is a long-standing problem in 
. statistical physics. Many of its conceptual aspects are similar to those encountered 
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in classical physics, but quantum mechanics poses the additional, but closely related, 
problem of decoherence of the quantum state of the system. As in classical physics, 
one needs to distinguish two easily confused cases: (1) The loss of information about 
5o i the state of a system due to its interactions with its environment, and (2) the loss of 

■ practically obtainable information due to the increasing complexity of its quantum 

^ , state. The first case is by now well understood through the separation of slow and 

fast degrees of freedom.^)"''') The growth of entropy of the subsystem (often involving 
only slow degrees of freedom and denoted by S) occurs because its state becomes 
5^ , entangled with the state |!P'e) of its environment, and the projection of the full density 

matrix p onto the system S implies a loss of information described by the "relevant" 
entropy 

Sroi = Tr [/9s In /9s] with ps = Tve p- (1-1) 

The second case, entropy growth in an isolated quantum system, is less well 
understood. However, it is an important issue both in cosmology, where the early 
universe is thought to have made the transition from a vacuum state to a thermalized 
state at the end of cosmic inflation, and in nuclear physics, where thermal matter is 
thought to be formed in collisions between two heavy nuclei moving at relativistic 
energies, both of them being in their quantum mechanical ground state. How can 
the evolution from a pure quantum state to a thermal ensemble be reconciled with 
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the unitarity of the S matrix, in other words, the tendency of, e.g., the Schrodinger 
equation to preserve the purity of the initial quantum state, and how can entropy 
growth be described in the absence of information loss to the environment? In classi- 
cal physics, the rate of information loss of nonlinear dynamical systems is described 
by their Kolmogorov-Sinai (KS) entropy, which is given by the sum over all positive 
Lyapunov exponents > of the system: 

SKS = J2^kO{Xk). (1-2) 

k 

The conditions under which the KS-entropy describes the growth rate of the entropy 
of a classical system have been widely explored (see e. g.^^)- There is no obvious gen- 
eralization of the concept of the KS-entropy to quantum systems, whose evolution 
is governed by a linear equation like the Schrodinger equation. Nevertheless, pure 
quantum states are known to evolve under their own intrinsic dynamics into state 
which can be, for most intents and purposes, be described as thermal ensembles. 
Compound nuclear states are a well known example. The theory of quantum chaos^^ 
addresses the question under which conditions highly excited eigenstates of a quan- 
tum system can be approximately represented as members of a thermal ensemble. ^'^^ 

Here we are addressing a more practical question: What is the time scale on 
which the state of a quantum system evolves from a simple structure easily recogniz- 
able as a pure state to a complex structure which is recognized by a typical observer, 
who himself is limited by the uncertainty principle, as an incoherent ensemble? Can 
the associated growth of (apparent) entropy be defined independent of the specific 
details of the measurement process and how can it be calculated? The starting point 
of our investigation is the expectation that, like in classical mechanics, the growth of 
entropy is governed by the dynamics of unstable modes, which amplifies uncertainties 
in the initial conditions. 

One may be tempted to think that the Wigner function, which furnishes a phase- 
space description of quantum dynamics, provides a useful starting point to address 
our question. However, being not positive definite, the Wigner function does not 
permit a probability interpretation. Also, the Wigner function is normalized such 
that the occupied phase space volume essentially stays constant with time. This 
shows up in the fact that expanding modes always are accompanied by contracting 
modes (classically, the positive and negative Lyapunov exponents for a conservative 
system come in pairs), see Fig.l for an example. The Husimi transform of the Wigner 
function, on the other hand, is positive semi-definite and admits a probability 
interpretation, but it also incorporates the limited ability of a typical observer to 
make measurements on the evolving system which can result in an increase of phase 
space volume, see Fig. 2. It thus forms an appropriate basis for our investigation. 

There are many discussions of entropy production in the literature which follow 
similar lines^^^"^^-* and address the same question for general or specific nonlinear 
dynamical systems. For example, Zurek and Paz^^^ use the example of the inverted 
(unstable) harmonic oscillator to show that the presence of a dissipative interaction 
with the environment leads to decoherence and the growth of the (von Neumann) 
entropy of a quantum system. Here we are mainly interested in applications to 
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problems in nuclear and particle physics, where interactions with the environment 
are either weak or absent, and the growth of entropy apparent to an observer is due to 
the rapidly growing internal complexity of the state of an essentially closed quantum 
system. Thus, our approach to the problem of entropy growth is different and does 
not involve interactions with an environment. As far as we know, this approach as 
well as its application to problems in particle physics has not been described before 
in the literature. 

In Sections [2] and [3] we investigate the Husimi function for two analytically 
tractable cases: the decay of an unstable quantum state and a standard toy model 
of reheating (or rather a variant called pre-heating) after cosmic inflation. We show 
that the growth rate of the entropy associated with the Husimi function approaches 
the classical Kolmogorov-Sina'i entropy in both cases. In Section 2] we study entropy 
production during the reheating (pre-heating) phase of cosmic inflation. In Section 
Owe compare the present entropy with other definitions such as the von Neumann 
entropy or the kinetic entropy. Finally, in Section El we discuss various possible 
applications of our formalism, including the production of entropy in the early stages 
of a relativistic heavy ion collision. 

§2. Entropy growth rate 

Following Zurek and Paz^^^ we choose the inverted harmonic oscillator as a 
simple example of an unstable mode in a quantum system: 

n = \f- \x'i'. (2-1) 

Here and below we distinguish quantum mechanical operators, such as the momen- 
tum operator p from classical quantities, such as the momentum p, by the caret 
symbol. We assume that the initial state of the system is given by a Gaussian wave 
packet of width sjh/uj: 



The Wigner function associated with the density matrix p is defined as 

W{p,x-t) = Jduei^P^ix-^imi^ + l)- (2-3) 

It contains the full quantum mechanical information of the system and provides 
a phase-space picture in accordance with the uncertainty principle. The Wigner 
function is easily seen to satisfy the normalization condition 



/ 



27rh 

Similarly, one finds that 



'^^'^''W{p,x;t) = Tt[p] = 1. (2-4) 



/ ^mP,^;t)? = Tr[p']<l. (2-5) 
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There is no simple relationship between higher moments of the Wigner function and 
Tr [/5"] for n > 2 . The Wigner function does not, in general, have a probabilistic in- 
terpretation, because it can take negative values as one easily confirms by calculating 
the Wigner function for an excited state of a harmonic oscillator. 

The time evolution of W{p, x; t) is determined by the equation of motion for the 
density matrix p: 

ih^^m = [^.m]- (2-6) 

For the initial pure state, (j2-2p the density matrix is /5(0) = |^o)(V'o|- As long as 
the Hamiltonian is at most quadratic in x and p or in the lowest order in /i, the 
Wigner transform of Eq. (|2-6p is equivalent to the Vlasov equation for W , dW/dt = 
{WjVFjp.B., where 7i{p,x) is the classical Hamiltonian. Then the Wigner function 
is constant along the classical path, 



xq cosh Xt + po/X sinh At, p = Xxq sinh At + po cosh At, 



(2-7) 



where {xq,pq) are the phase space variables at t = 0. The time-dependent Wigner 
function is then easily obtained as. 



W{p, x;t) = 2 exp 



1 



2 exp 



Kip,x;t) 

h 



(2-8) 



with 



K{p,x;t) 



P2_ 
X 



{acosh2Xt + 6) + Xx'^{a cosh2Xt - 5) - 2apxsmh2Xt. (2-9) 



The parameters a and 6 are defined as 

A2+a;2 



> 1, 



A^ 



2Xoj ~ 2Xoj 
The Husimi function is defined as the Gaussian smeared Wigner function: 



(2-10) 



HAiP,x]t) 



dp' dx' 



■ exp 



exp 



hA{t) 



^{x-x'f 
n 

K{p,x;t) + ^ + Ax^ 



W{p',x';t) 



where 



A{t) = 2(cjpcosh2At + 1 + 65') 
and the parameters p and 5' are defined as 



P 



Z\2 + A2 



> 1 



5' 



2Z\A - 2Z\A 
In our convention the Husimi function is normalized according to 



dpdx 
27rh 



W{p, x; t) 



dpdx 

2TTh 



HA{p,x;t) = 1. 



(2.11) 
(2-12) 

(2-13) 
(2-14) 



Towards a Theory of Entropy Production in the Little and Big Bang 



5 



In contrast to the Wigner function, the Husimi function is always non-negative, 
because it can be expressed as the expectation value of the density matrix in a 
coherent state: 

HA{p,x;t) = {zA\m\zA), (2-15) 



where a/\|z^) = za\za) with 



aA = — , (2-16) 

V2hA 




Fig. 1. The Wigner function (|2-8[) for the unstable osciUator ai t = and t = 2/X for a — 1. The 
horizontal axis denotes the scaled position ujx; the vertical axis represents the momentum p. 




Fig. 2. Husimi function (|2-11[) for the unstable oscillator at t = and t = 2/A for p = a = 1. Note 
that the extent of the distribution in the off-diagonal direction (p — ujx) does not shrink beyond 
the resolution limit set by the Gaussian smearing introduced by the Husimi transform. 
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Because it is non-negative, the Husimi function can be used to define a coarse 
grained entropy of the quantum state, first introduced by Wehrh^^^ 



To simphfy notation we introduce 



dpdx 

27Th 



Ha{p, x; t) In Ha{p, x; t) 



(2-17) 



L{x,p,t) = K{x,p,t) + ^ + Ax'^ 



(2-18) 



and get 



Sn,Ait) = - 





dhj ^/A(i 
dp dx 



2-K 



HA{p,x;t) 



In 



+ 1 



- In — ^ + 1 
2 4 



(2-19) 



The Wehrl entropy is a measure of the complexity of the state of the system, as one 
can see as follows. The volume of support in (2i5)-dimensional phase space of the 
Wigner function of a D-dimensional system in a pure quantum state is always equal 
to = (27r/i)^, independent of time. This property is mathematically expressed 
by the fact that the square of the Wigner function for any pure quantum state is 
equal to one: 



1. 



(2-20) 



This property is a consequence of the relation Tr [fp] = Tr [p] = 1 valid for the density 
matrix of a pure quantum state. For an unstable or chaotic system, the shape of the 
volume in which VF(p, x; t) is significantly different from zero becomes more and more 
elongated and irregular as time progresses. If we divide phase space into a regular 
grid of cells of volume , the Wigner function will take on substantially nonzero 
values in an increasing number of phase space cells. If we define the complexity of 
the state of the system as the number of cells in which the Wigner function averaged 
over the cell takes values larger than a given lower bound, the complexity will grow 
with time as the support of the Wigner functions becomes more and more irregular. 

The Husimi function measures this growth in complexity by smearing VF(p, x; t) 
locally with a minimum uncertainty Gaussian of width . The volume of support 
of the Husimi function in phase space is thus equal to the number of phase space cells 
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M,^ I -J^HT-\p,x-,tf<l. (2-21) 



"touched" by the Wigner function in a significant way. The growth in the volume 
of support of the Husimi function in phase space can be expressed in terms of its 
second moment 

2'Kh 

The inverse of M2 has been proposed as a convenient and easily calculable measure 
of the complexity of the state of a quantum system. ^^^'^^^ Its great disadvantage, for 
our purposes, is that it is not a measure of the coarse grained entropy of a quantum 
system and does not approach the von Neumann entropy for a system described by a 
thermal ensemble. The Wehrl entropy, 5'h,A) on the other hand, shares in the ability 
of the Husimi function to trace the complexity of a quantum state, and it permits an 
interpretation as the entropy of a quantum system. We thus focus our investigation 
on it. 

The rate of growth of the Wehrl entropy 

dS}i A f dpdx dH/\ . d f dpdx f dp dx dH^ , 

lnii/\ + — / -——Ha= / -r-^: ^—IiiHa 



dt J 2-Kh dt dt J 2'Kh J 2'Kh dt 

^ \ap sinh2At t->op ^ ^2.22) 
ap cosh 2\t + l + 55' ' ^ ' ' 

asymptotically tends to the growth rate of the unstable mode, i.e. to the positive 
Lyapunov exponent of the classical Hamiltonian. This result has several noteworthy 
features. First, the entropy growth rate tends to the classical Lyapunov exponent 
very rapidly, on the time scale of itself. Second, the value of the growth rate for 
t ^ A~^ is independent of the smearing parameter A. Although the absolute value 
of the Wehrl entropy ()2-19p depends logarithmically on A via p and 5' , its growth 
rate does not. 

Let us repeat and expand a bit the basic physics idea behind the analysis just 
presented: Any kind of measurement results in some coarse graining, which in turn 
can result in entropy production for a quantum system. The question is whether 
one can define entropy independently of the specifics of a measurement, only taking 
into account the principal limitations imposed on any measurement by quantum 
mechanics. The Husimi function provides an answer, as it incorporates the fact that 
any measurement (where here we only regard those which can be parametrized by 
a width A) must fulfill the uncertainty principle. In general, the value of the coarse 
grained entropy will depend on the value of A and thus will be ambiguous. The 
central result obtained here is that the entropy growth rate at sufficiently long, but 
not too long, times is independent of A. We conjecture, but have not proved, that 
the entropy growth rate is independent of other parameters one could introduce to 
specify an actual measurement and that the quantity dSn^^/dt is a well-defined, 
physical property of the system and not of the measurement process. We note that 
the connection between the exponential rate of growth of the phase space volume 
occupied by the Husimi function and the classical Lyapunov exponents was pointed 
out before, e. g., by Toda and Ikeda.^^^ These authors, however, did not draw 
the connection between the growth rate of the Wehrl entropy and the Lyapunov 
exponents of the classical theory. 
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Next we consider a system with a large (possibly infinite) number of unsta- 
ble modes. Assuming no interactions among these modes, the Hamiltonian can be 
written as 

^ = E^(^5fc-^M) (2-23) 



2 

k 



Assuming that the initial state is a product of Gaussians of the form of Eq. (j2-2p the 
Wigner function and the Husimi function of the complete system are given by 

W{{pk},{xk};t) = l[w^''\pk,xk;t) (2-24) 

k 

HAi{Pk},{^k};t)=l[H^^\pk,Xk;t) (2-25) 

k 

Using the normalization condition for the Husimi functions of the individual modes 

it 

5H,A(t) = (2-26) 



(k) 

it is now easy to show that 



k 

The growth rate of the total Wehrl entropy is then 



c^'S'h.a _ ■sp Afc sinh2Afct 

dt ~ ^ cosh2Afct+ (1 + 

k 

which, for long times, tends to the Kolmogorov-Sinai entropy of the system: 

dSii^A t^oo 



dt ^ 



^Afc. (2-28) 



The limit is approached in such a way that the large contributions (large Afc) are 
reached rapidly while the small ones take longer. The dependence on the details of 
the Gaussian smearing disappears exponentially with time. 



§3. "Roll-over" transition 



We next consider a system with infinitely many degrees of freedom, namely, 
relativistic quantum fields with dynamically unstable modes. Here we will see that, 
while the entropy grows linearly with time for unstable modes, it oscillates around 
a constant for stable modes. We will also find that the entropy growth for a trans- 
lationally invariant instability rate is an extensive quantity, i. e. it grows with the 
volume. 

For the purpose of analytical tractability we choose the scalar field in (1+1) 
dimensions with the tachyonic Lagrangian (we set ^ = 1 in this section.) 
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The modes with momentum \p\ < // are unstable, i. e. their amphtudes grow ex- 
ponentially with time; the modes with \p\ > fi are stable and exhibit oscillatory 
behavior. The Wigner functional W[n{x)^<P[x)] for the field ^ and its canonical 
momentum 77 = d<P/dt can be defined in complete analogy to the Wigner function 
of a quantum system with many degrees of freedom^^) as: 



W[n{x),^{x)]t] = / Vip{x) e-»/'^^^(^M^) 



1 



1 



x{${x) + -((C(x)l f){t) \<^{x) - -ip{x)) 



(3-2) 



The Wigner functional for t > can be represented in terms of the Fourier 
modes 77p, ^p. We note that the Wigner functional is a constant along the classical 
trajectory (^p(i), 77p(i)): 



W^[77,<7>;t] = VF[{770},{^>0},t = 0] 



(3-3) 



where and 77^ are the initial Fourier components. If we assume that at i = the 
scalar field is in the vacuum state corresponding to mass m, this property allows to 
write the Wigner functional as: 



W[n,^]i\ = C exp 
where for < p < //: 

<Pp = ^p{t) cosh Apt 
rO 



dp 
2^ 



|770|2 ^ 



(3-4) 



77p(t) 



sinh Apt 



77p = 77p(t) cosh Apt - Ap (Pp(t) sinh Apt 



with Ap = — p'^, and for p > /x: 



(Pp = (Pp(t) cos Wpt 



77p(t) 



sin LVpt 



77° 

p 



Upit) cos ujpt + u}p 'Pp{t) sin ujpt 



(3-5) 
(3-6) 

(3-7) 

(3-8) 
(3-9) 



with Ep = Y^p2 + and 



Ldr, 



The evaluation of the Wigner and Husimi functions proceeds in complete analogy 
with the case of a single unstable mode discussed in the previous section. The final 
result is: 



HA[n,'P;t] = H 



exp 



R{np,^p-t 

Ap{t) 



n 



ip(t) 



exp 



R{np,^p-t) 



Ap{t) 



(3-10) 
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Fig. 3. Contribution of an unstable mode, shown as solid line to the Wehrl entropy compared with 
the contribution of a stable mode, shown as dashed line. The contribution of the unstable mode 
grows linearly with time, while the contribution of the stable mode oscillates around a constant 
value. The curves are for the parameters \p = ujp = 1 and Ep = A = 2. 



with 



\n, 



R{np,^p]t) = — -^^ — (dp cosh 2Apt + (5p) + Ap|<?p|^((7p cosh 2Apt — (5p 



A 



-Op{nl^p + Hp^l) sinh 2Apt + 



A 



p\ 1 



R{np, (pp] t) = — - — (cjp + 6p COS 2u!pt) + u!p\(pp\'^ {dp — 6p cos 2u!pt) 



177 P 

+6p{n;^p + 77p^;) sm 2iOpt + + A\^p\^; 



(3-11) 



with 



Apit) 
Ap{t) 



— Or, cosh 2Ar,t + 2 + 6r, 



XpA 



dp + 2 + 6p 

ujpA 



p 

. A^ 



cos 2ijjpt 

UJnA 



an 



(Tr, 



2XpEp 



> 1, 



> 1, 



_ Ap - £/p 

Op 



2XpEp 



5n 



2ujpEp 



Like in the single mode case, the Wehrl entropy is given by 

DnD<p 



2tt 



(3-12) 



(3-13) 
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V 



dp 
2^ 



InMl + l 



+ V 



dp 
2^ 



2 4 



(3-14) 



where V is the quantization volume. In analogy to (|2-22p . the growth rate of the 
entropy is given by 



dt 



V 



dpcjp(Z\2 + A2)sinh2Apt 



V 



27r 
f dp 



Ap{t)A 



+ V 



dphi^l 



zl^) sin 2u)pt 



bl>M 



27r 



Ap{t)A 



(3-15) 



It is instructive to compare the integrands in (j3-14p for a typical stable and 
unstable mode. Figure [3] shows such a comparison. One clearly sees that the contri- 
bution to 5h,a of the unstable mode (shown as solid line) grows linearly with time, 
while the contribution of the stable mode (shown as dashed curve) does not grow 
and only oscillates slightly around its initial value. It is also noteworthy that the rate 
of entropy growth (|3-15p is proportional to the volume V , implying that the entropy 
growth rate is an extensive quantity. This means that the evolution of the scalar 
field is characterized by a finite growth rate of the entropy density S/V . Generally, 
the entropy growth rate will be an extensive quantity, if all modes below a certain 
momentum scale are dynamically unstable. 



§4. Entropy production in the "Big Bang" 



A different case where the dynamics of the system is thought to lead to coarse 
grained entropy production is the "(p)re-heating" after inflation problem. ^'^^ In this 
case one often considers the following toy model: A scalar field x (representing 
matter degrees of freedom) interacts with the inflaton field ^ via the Lagrangian 

Because the infiaton field has the same value everywhere in space, the problem is 
translationally invariant and can be separated into independent modes by using the 
Fourier expansion (we set h = 1 m. this section): 

where R{t) is the scale factor describing the expansion of the universe. After infiation, 
the inflaton field ^ oscillates with a large amplitude and can be described, to good 
approximation, as a classical field: 

^(t) w <^o cos(a;i) (4-3) 

After rescaling the physical time variable as r = ytjtc with a characteristic cos- 
mological time constant tc, the equation for the quantum operators describing the 
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individual Fourier modes of the matter field then takes the form of a Mathieu equa- 
tion (see e.g}^^): 

+ f k2 + 1^ cos(2..t)^ Mr) = 0, (4-4) 



where A is the self-interaction constant of the inflaton field and = fc^/(A<?g). The 
solutions of this equation are the Mathieu sine and cosine functions S{a, q; ojt) and 
C{a,q;u!T) with a = and q = —g^/A\. We will henceforth drop the parameters a 
and q. Asymptotically, the Mathieu functions can be represented as 

C{ujt) m e''^ cos{ujT + ac{T)) , S{u)t) ^ e''^ cos(a;r + as(r)) (4-5) 

with /L( > and ac{0) = 0, as{0) = — 7r/2. For large times, the constancy of the 
Wronskian {CS — SC)/uj = 1 ensures that the phases of both solutions approach 
each other with exponential precision, according to 



\a,{T) - as{r)\ e-'^\ (4-6) 

In order to explore the consequences of this contraction in phase angle, we 
choose the phase-amplitude representation of phase space for each mode: (X^, P^) = 
(rfc cos Ofc, wrfc sin Qffc). Thus we have in this case as conjugate variables the squared 
amplitude r^/2 and the phase ak- Because in our case, both Xfc(r) and Pkir) asymp- 
totically oscillate with the same phase adr) « ces{T), the uncertainty of the particle 
number (see Eq. (238) in^o)) 

is equal to rather than \/Nk- Thus we have 

iANk)\Aakf = NliAauf > \. (4-8) 

Although it is quite commonly used, the status of ()4-8p is not unproblematic (see, 
e.g., the discussion in^^^). In general, the expectation value of the commutator 
[iVfc,Q!/c] is state-dependent. Our states (the solutions of the Mathieu equation) 
are characterized completely by A^^ and and for them Eq. (|4-8p holds. Thus we 
here regard Nk and at as canonical conjugate for the states we consider. As a 
result Eq. (j4-8p generates coarse-graining, which in turn can be linked to entropy 
production. Therefore, the Wehrl entropy defined below is a sensible way to define 
entropy for our setting. 

We also note that the 27r-periodicity of implies that r|,/2 is an integer. How- 
ever, for large enough r, this quantization becomes irrelevant, and can be treated 
as a continuous variable. 

In the following we drop the mode index k. We assume that the wave function 
!P'(a, r) is nearly independent of the phase angle initially and then contracts according 
to (l¥6D : 

^{a,T) ^ (fr(r)V^)-^/%xp (-^^) , (4-9) 
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where Tr{T) = vr exp(— 2^r). We shall drop the argument of vf in the following, but we 
assume r to be large enough so that vf (r) <C n. It is now a straightforward exercise 
to calculate the Wigner function: 



W{a, n, t) 



j d/3e*"/3g^p/ ^\ ^2e-" * , (4-10) 

where n = is an integer. For continuous variables r would be e^^. The Husimi 
function is similarly obtained as 

HMa,n,r) = £ ^ J] e-^("-"')^-^ W^(a', n', r) 



2V^^ f o?A + n^vi-^ 



...... (4-11) 





Fig. 4. Wigner function (|4-10[1 for the inflationary reheating model at [ir — 0.5, 1.5. The horizontal 
axis denotes the phase angle a G [— 7r,7r]; the vertical axis represents the squared amplitude 
e [0,10]. 



Finally, we obtain the Wehrl entropy: 
Sr,a(t) = - 



2fiT + const. 



In — , ^ ^ + - 



(4-12) 



Once again, the rate at which the Wehrl entropy increases is given by the classical 
Lyapunov exponent 2/i. Taking several independent modes into account one obtains 
as before the sum of all positive Lyapunov exponents, i.e. the Kolmogorov-Sinai 
entropy. 
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Fig. 5. Husimi function (|4-lip for the inflationary reheating model at fir = 0.5, 1.5 and for A — 2. 
The horizontal and vertical axes are chosen as in Fig. |4l 

§5. Relation to Other Definitions of Entropy 

5.1. Thermal equilibrium 

We next discuss the relation between the Husimi- Wehrl entropy discussed in the 
previous section and the von Neumann entropy in thermal equilibrium. The simplest 
system to consider is a harmonic oscillator, 

W = ^ + ia;V . (5-1) 

It is also a quite generally applicable, because a quantum field at weak coupling or low 
excitation can be considered as an ensemble of infinitely many harmonic oscillators. 
The occupation probability of the eigenstate |n) is given as, 

wn = e-^P'^/Zp , (5-2) 

oo 

Zf3 = Y^ e""^^ = (1 - e"^'^'")-^ . (5-3) 

n=0 

where P = 1/T denotes the inverse temperature, and is the partition function. The 
eigenstate |n) corresponds to the n-th excited state in quantum mechanics and to 
the n-particle state in a Fourier mode u = Ep in the case of a field theory. The von 
Neumann entropy is calculated to be 



oo 



vN 



n=0 



w„ In w„ 



f3huj 



-n In n + (n + 1) ln(n + 1) . 



ln(l 



-i3nuj\ 



(5-4) 
(5-5) 
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In the second line, the average number of excited quanta (or the average number of 
particles in the field theory) is used, 

^ n=0 

The thermal Wigner function is given by^^^ 

W{z) = Bpexp{-Bf3Zz) , (5-7) 

where = 2 tanh(/3^/2) = l/(n+l/2), and we have applied a complex representa- 
tion of {x,p) as z = {ujx+ip) / V2hw . Although the Wigner function is, in general, not 
positive definite, it is positive around the thermal equilibrium and can thus be used 
to define an entropy, which we here call the Wigner- Wehrl entropy S'w ■ By using the 
phase space measure, dP = dxdp/27rh = dzdz/TT{= (i(Rez)(i(Im2;)/7r) = d{zz)d6 /2tt, 
the Wigner- Wehrl entropy is obtained as, 



Sw = -j dPW{z)\nWiz) = - j ^^^Bpe-^'p'^nniBpe 



= l-\nBp = l + \ii{n+^ , (5-8) 

where denotes the argument of z. 

The Husimi function is the expectation value of the density matrix with the 
coherent state defined as, 

oo 



n 

exp(zat) |0) = e-^"^/2 ^ |n) , (5-9) 



where we have chosen the natural width Z\ = a; in Eq. (|2-16p for simplicity. The 
Husimi function in thermal equilibrium is given as. 



H{z) = {z\ Ah k) = ^ E 

n=0 

1 

^ exp 



Z(3 ^ n\ 



-zz 1 



{l-e-f'^)]=A(,e^p{-Apzz), (5-10) 



where pth = SJ^Lo W) '^n {n\ denotes the density matrix in thermal equilibrium, 
and Ap = 1 — e~^^ = l/(n -|- 1). The Husimi- Wehrl entropy then is obtained as 

Sn = - J dPH{z) \nH{z) = - J ^^^^Ape-^^'' \n{Ape-^0'') 

= l-lnAf3 = l + ln{n+l) . (5-11) 



The Wigner- Wehrl and Husimi- Wehrl entropies, 5w and Su in Eqs. (j5-8p and 
(j5-lip . respectively, differ from the von Neumann entropy, 5vN in Eq- (!5-4p . at low 
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temperatures (T < hio), while all three entropy definitions approach the same asymp- 
totic value, 1 + Inn ~ 1 — In f3huj, at high temperatures (T ^ huj). As discussed 
in Ref. 24), the Wehrl entropies are always greater than or equal to the von Neu- 
mann entropy (Wehrl-Lieb theorem), and in the present case, we have shown that 
the difference is at most one unit as far as we adopt the natural width. The tem- 
perature dependence of the von Neumann entropy and the Husimi- Wehrl entropy 
for the harmonic oscillator are plotted in Fig. [6l Since the average occupation num- 
ber of a quantum field mode in thermal equilibrium grows with temperature as 
fi T/huj ^ 1 (for T » hw) we can surmise that, in general, Su will be a good 
approximation to the physical entropy for field modes that are highly occupied. In 
the application to the equilibration phase of relativistic heavy ion collisions, which 
we will discuss in the next section, this condition is satisfied, because in the domain 
of interest, the so-called "glasma" phase, the average occupation number of field 
mode is ?!. 1/q<j ~ 3, where the difference between Sh and S^N is less than 10 
percent. 



4 




The difference between the Husimi- Wehrl entropy and the von Neumann entropy 
is caused by the energy fluctuation in the coherent state. In order to demonstrate 
this point, we consider the case where the partition function is given in the coherent 
state representation as follows, ^^-^ 



2/3 = J dr{z\ e-^^ \z) = J drexp - dp'Tip^iz) 



d(3 



\n{z\ e-^^ \z) = {z\ He-^'^ \z) / {z\ e"^^ \z) 



{z\n\z)-PaUz) + 0{f3^) 



(5-12) 



(5-13) 
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where 7ip{z) is the thermally distorted expectation value of the Hamiltonian, and 



(y^{z) represents the energy variance <t^{z) = {z\ 7i 



{z\7i \z) . From general 



thermodynamic relations, the (von Neumann) entropy is calculated as 

d 



OF d 

{z\ e-/^« 



dr- 



[(3Hp{z)+\TiZp\ . 



(5-14) 



The difference from the Husimi-Wehrl entropy is evaluated to be, 



Sn-S = j dr H{z) 
~ j drH{z) 



(5-15) 



This relation shows that Husimi-Wehrl entropy agrees with the von Neumann en- 
tropy when the energy fluctuation is much smaller than the temperature. Systematic 
corrections are possible when we can evaluate the expectation value of Ti"". Specif- 
ically, in the case of the harmonic oscillator, the evolution of fi^z) with /? and the 
energy variance are known as, Ti-piz) = {z\ H \z) exp(— and cr'^{z) = {hio)'^zz, 
allowing us to obtain explicit correction terms which exactly reproduce the von Neu- 
mann entropy, if this is desired. 

5.2. Dynamical versus kinetic regime 

For systems which permit a quasi-particle description, the Kadanoff-Baym for- 
malism^^) can be used to derive the kinetic theory from quantum field theory. The 
derivation usually starts from the single-particle Wigner function 



Wsp{p,x;t) = / dueftP"(^(x 



(5-16) 



and performs the Moyal (gradient) expansion to define the phase space distribution 
f{x,p,t). One can then define the kinetic entropy via the standard phase space 
integral (see e. g.^^^): 



S-kh 



|f|f[-/ln/ + (l + /)ln(l + /)]. 



(5-17) 



(For fermions, the plus signs are replaced by minus signs in the second term in square 
brackets.) The kinetic entropy so defined agrees with the von Neumann entropy in 
equilibrium and satisfies the i?-theorem for processes relaxing to or close to the 
thermal equilibrium state. An extension to include memory effects in the kinetic 
equations is also possible. The Kadanoff-Baym formalism was recently applied 
by Nishiyama to study entropy production in the kinetic regime of the scalar <P^- 
theory.^^) Its extension to nonabelian gauge fields can be found in the review of Elze 
and Heinz;^^) the application to relativistic nuclear collisions was explored by Gelis 
et a/.30) 
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The Kadanoff-Baym formalism makes various approximations, some of which 
break time-reversal invariance, in order to derive a kinetic equation: (i) A truncation 
of the Feynman diagrams at the two-particle irreducible (2PI) level is imposed. This 
truncation implies an information loss, because two- and more-particle correlations 
in the quantum state are neglected, (ii) Since it is impossible to include all the 
diagrams of the 2PI type, one must further make a skelton expansion and take only 
diagrams with a limited number of loops, {iii) Finally, the gradient expansion is 
made to break the space-time nonlocality of the quantum theory. This amounts to 
a Markovian approximation. 

Here, on the other hand, we are interested in the passage from the dynamical 
(Liouville) regime of the quantum field theory to the kinetic regime. This problem 
arises when the system of interest starts in a highly excited, but still nearly coherent 
quantum state and then evolves toward equilibrium. In such situations, one cannot 
neglect the many-particle correlations which are encoded in coherent field configu- 
rations, because it is precisely the loss of coherence of such configurations due to 
their own nonlinear dynamics, which leads to the initial growth of entropy. The 
Kadanoff-Baym approach to kinetic theory is not applicable to this domain. 

It is not clear at all, whether a unique definition of the concept of entropy can 
be given, which permits the description of the transition from pure quantum state to 
thermal equilibrium and applies to all possible systems. The definition of entropy is 
firmly linked to the concept of coarse graining of the quantum state and thus depends 
on the specific chosen procedure. The Husimi distribution is one possible choice; it 
is motivated by the idea that coherent field configurations play an important role 
during the transition from the (nearly) pure quantum state to the kinetic regime. We 
emphasize, however, that we do not assign special importance to the absolute value 
of the Husimi- Wehrl entropy — since its precise value depends on the arbitrary 
smearing parameter A, this would be difficult to justify — but only to its growth rate 
dSn/dt, which is independent of A in the long-time limit. It would be interesting to 
study how the growth rate of the Husimi- Wehrl entropy compares with the growth 
rate of the kinetic entropy, as a system enters the kinetic regime. We leave this 
question to a future study. 

§6. Relativistic heavy ion collisions 

Much experimental evidence exists showing that a thermal quark-gluon plasma 
is formed in collisions between two heavy nuclei at center-of-mass energies on the 
order of 100 GeV per nucleon.'^^) The total entropy per unit of rapidity, dS/dy, 
of the final state in such a reaction of two Au nuclei has been estimated to be 
approximately 5,000.^^^'^^^ How, and how rapidly, is this entropy produced? Exper- 
imental observations, in combination with hydrodynamics simulations, suggest that 
the thermalization time is very short, of order 1 fm/c or less.^^^ A fraction of the 
final entropy can be produced simply by decoherence of the initial quantum states 
of the colliding nuclei, on a time scale Tdec ~ l/Qs; where Qs is the saturation scale 
of the nuclear parton distributions. ^^^'^^^ 

The remainder of the entropy can be generated at various stages of the heavy-ion 
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collisions: One source which is certainly relevant are nonlinear interactions among 
the gluon fields liberated in the decoherence process, but the precise mechanism 
responsible for the entropy growth and thermalization is still unclear. Suggested 
processes include the Boltzmann cascades of partons,^^^"^^) the nonlinear dynamics 
of classical color fields, ^'^^"^^^ the decay of unstable chromo-magnetic fields, ^^^'^^^ and 
plasma instabilities in the longitudinally expanding matter. '^^^"^^^ Another relevant 
source might be the viscosity terms in viscous hydrodynamics. ^^^"^^^ The hadronic 
phase and processes like jet fragmentation^^) are promising sources. Here we con- 
centrate, however, on the initial phase of the heavy-ion collision. 

All approaches cited above are based on classical or quasi-classical pictures, 
in which the production of entropy is mostly assumed (as in parton cascades) or 
based on plausibility arguments (as in simulations of classical Yang-Mills fields), 
but not rigorously calculated. As far as we know, no comprehensive formalism 
has ever been proposed which would, even in principle, describe the calculation of 
entropy production starting from the initial quantum states of the colliding nuclei. 
Fukushima et al.^^^ analyzed the quantum fiuctuations around classical glue fields 
in the earliest phase of the collision, but did not address the problem of entropy 
production. Their analysis is also problematic in our context, because it assumes 
that the collision time is exactly zero. 

In principle, the Husimi function approach described in the previous sections 
provides the missing link in the formulation of a comprehensive approach to entropy 
production in nuclear collisions. The lattice regularized Yang-Mills equations de- 
scribing the dynamics of classical color fields were found to be strongly chaotic, ^^^^^''^ 
and the KS-entropy of the classical Yang- Mills field was shown to be a thermodynam- 
ically extensive quantity.^^^^ The general formulation of the dynamics of quantum 
fields based on the Wigner functional has been developed by Mrowczynski and one 
of the present authors. An extension of the classical lattice Yang-Mills theory to 
include Gaussian quantum fluctuations around the classical link variables was pro- 
posed by Gong et al.^'^^ Based on these developments one can construct a Husimi 
functional for the lattice Yang-Mills held, using either the Wigner functional or the 
Gaussian variational method, and thereby link the classical calculation of the KS- 
entropy density to the entropy density growth rate of the quantized gluon field. 

§7. Summary 

We have proposed a method to determine the rate of entropy production for 
quantum field theoretical systems like colliding nuclei or the early universe. This 
method employs the Husimi function to implement the amount of coarse graining 
required by the uncertainty principle. For the entropy we adopt the definition of 
Wehrl. We then flnd that for large enough times the rate of entropy growth is given 
by the analog of the classical Kolmogorov-Sinai entropy. This observation agrees with 
results from various discussions found in the literature. We illustrated these general 
ideas for a few concrete examples. In Sections [2] and [3] we analyzed the situation for 
a scalar field in an inverted harmonic oscillator potential. The standard toy-model 
for pre-heating after cosmological inflation was analyzed in Section [H In Section [5] 
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we compared the Wigner-Wehrl and Husimi-Wehrl entropy with other definitions of 
entropy, such as the von Neumann entropy and kinetic entropy. 

Based on these observations and examples we have shown that it is possible 
to define the conceptual framework for a calculation of the entropy growth rate in 
heavy ion collisions which rests on a much firmer theoretical basis than previously. 
In Section [6] we outlined the steps necessary to carry out such a calculation. Its 
realization will require a prolonged and extensive effort. The initial conditions for 
the gluon Wigner distribution must be derived from the gluon wave functions of 
the colliding nuclei, which is only known approximately within certain models (see 
e. g.)^^^). The time evolution of the Wigner and Husimi functions for the Yang- 
Mills field will be strikingly different from that of the simple cases analyzed here, 
in which different modes did not couple. The time evolution as well as the Husimi 
transform will have to be done numerically, due to the complexity of the problem. 
The accuracy and precision of the applied numerical methods will have to be tested 
carefully. We should confirm, e.g., that the growth rate of the entropy density 
becomes independent of the number of lattice points used to discretize space for the 
quantized Yang-Mills field. However, all of these problems are technical rather than 
conceptual and will hopefully be resolvable in due course. 
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